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Abstract 

We show how to evaluate tensor one-loop integrals in momentum space avoiding 
the usual plague of Gram determinants. We do this by constructing combinations of n- 
and (n — l)-point scalar integrals that are finite in the limit of vanishing Gram deter- 
minant. These non-trivial combinations of dilogarithms, logarithms and constants are 
systematically obtained by either differentiating with respect to the external parame- 
ters - essentially yielding scalar integrals with Feynman parameters in the numerator 
- or by developing the scalar integral in Z) = 6 — 2e or higher dimensions. An ad- 
ditional advantage is that other spurious kinematic singularities are also controlled. 
As an explicit example, we develop the tensor integrals and associated scalar integral 
combinations for processes where the internal particles are massless and where up to 
five (four massless and one massive) external particles are involved. For more general 
processes, we present the equations needed for deriving the relevant combinations of 
scalar integrals. 



^Address after 1 January 1997, Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, 
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1 Introduction 



One of the most important ingredients in the search for "signals" of new phenomena in 
high energy particle physics experiments is a precise knowledge of the expectations from 
standard physics; the "background". Usually this involves perturbative calculations of dif- 
ferential cross sections within the standard model. Many such radiative corrections have 
been carried out and require the evaluation of one-loop integrals which arise directly from a 
Feynman diagrammatic approach. Often these integrals need to be performed in an arbitrary 
dimension in order to isolate any infrared and ultraviolet divergences that may be present 
The basic one-loop tensor integral in D dimensions for n external particles scattering 
with outgoing momenta pj, n internal propagators with masses Mj and m loop momenta in 
the numerator can be written, 



In 



where m = 1, . . . , n and, 



(£2 _ Mf){{i + q,y - M|) ■■■{{£ + - M^) ' 



i=i 

The scalar integral is denoted I^[l]. In the standard approach to such integrals one 
utilises the fact that the tensor structure must be carried by the external momenta or the 
metric tensor g'^'^. For example, the simplest non-trivial tensor integral in D dimensions has 
a single loop momentum i^. It reads, 

n-l 

i^n^E^p,", (1-1) 

using momentum conservation to eliminate one of the momenta. The formfactors Cj are 
determined by multiplying both sides by all possible momenta and rewriting i.pi as a 
difference of propagator factors, 

^■P^ = I [{{£ + q^f - Ml,) - {{i + - M^) + {Ml, -Mf + qU - q^) 

thereby reducing the tensor integral to a sum of scalar integrals, 

= I [ln-f'\l] - + (Ml, -M^ + qU - ql)m 

Here, I^1^J^[1] represents the "pinched" loop integral with the (n — 1) propagators remaining 
after the ith propagator factor has cancelled. 
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The formfactors are then obtained by algebraically solving the system of {n — 1) equations. 
This introduces the {n — 1) x [n — 1) Gram determinant, 

An = det{2pi -pj), 

(where i and j run over the (n — 1) independent momenta), into the denominator. Each 
formfactor is a sum over the scalar integrals present in the problem multiplied by a kinematic 
coefficient that may be singular at the boundary of phase space where the Gram determinant 
vanishes. Typically, 

Cj ~ + 1^ Prn , (1.2) 

where the sum is running over the m possible pinchings and where a and jSm are coefficients 
composed of the kinematic variables. Since, in many cases, the formfactors Cj are actually 
finite in this limit, there are large cancellations and there may be problems of numerical 
stability. 

The basic approach has been modified in a variety of ways, including the introduction of 
a system of {n — 1) reciprocal vectors ff (and the associated second rank tensor w^'^ playing 
the role of g'^'^) to carry the tensor structure 0, || where, 

,,/^ — ,Pl...p»-lMPi+l...Pnp /A 
~ ^ p^l 

SO that, 

Vi-Pj = Sij. 

This simplifies the identification of the formfactor coefficients, but does not eliminate the 
Gram determinants. In fact, in both approaches, the number of Gram determinants gener- 
ated is equal to the number of loop momenta in the numerator of the original integral. 

A different approach has been suggested by Davydychev ||^, who has identified the form- 
factors directly as loop integrals in differing numbers of dimensions and with the loop prop- 
agator factors raised to different powers. Tarasov has obtained recursion relations for 
one-loop integrals of this type, so that a complete reduction is possible. However, in relating 
the formfactor loop integrals to ordinary scalar loop integrals in 4 (or close to 4) dimensions, 
the Gram determinant once again appears in the denominator as in equation ( |1.2|) . 



Finally, Bern, Kosower and Dixon have used the Feynman parameter space formulation 
for loop integrals to derive explicit results for the scalar integrals including the scalar pen- 
tagon 0. The formfactors of the momentum space decomposition are directly related to 
Feynman parameter integrals with one or more Feynman parameters in the numerator. One 
can see this by introducing the auxiliary momentum V^, 



n-l 



= - 5: x.+igf , (1.3) 
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so that after integrating out the loop momentum, the tensor integral for a single loop mo- 
mentum in the numerator can be expressed in terms of the external momenta gf 



n-1 
i=l 

Here, represents the scalar integral with a single factor of Xi in the numerator. By 

comparing with equation (|1.1|), we see that, 

n 

Cj = — E i-^i]- 

i=j+l 

Differentiating with respect to the external kinematic variables, yields relations between 
integrals with polynomials of Feynman parameters in the numerator and the usual scalar 
integrals. Once again, the Gram determinant appears in the denominator, and the final result 
for the formfactor cj combines n-point integrals with the pinched [n — l)-point integrals as 
in equation (|1.2|) . 



The presence of the Gram determinant is, in some ways, no great surprise. In the 
limit An 0, the (n — 1) momenta no longer span an (n — l)-dimensional space, the 
(n — 1) equations of the Passarino-Veltman approach are no longer independent and the 
decomposition is invalid. Stuart []TU[ has made modifications to the basic approach to account 



for this, the main observation being that for A„ = 0, the scalar ra-point integral can be written 
as a sum of scalar {n — l)-point integrals. As a consequence, there are large cancellations 
between scalar integrals with differing numbers of external legs in the kinematic limit of 
vanishing Gram determinants. For loop corrections to processes such as quarkonium decay, 
where the two heavy quarks are considered to travel coUinearly and share the quarkonium 
momentum, one can eliminate the Gram determinant singularities completely using the 
method of Stuart pO[| . 

However, for more general scattering processes where the coUinear limit may be ap- 
proached, but is not exact, the numerical problems as A„ remain. In this paper, we 
wish to address the problem by combining the scalar integrals into functions that are well 
behaved in the A„ limit |]TT[, so that the formfactors are given by. 



9 - «^ ( I^[l] + E P'mln^T\l] + finite. 



A,, 



m=l 



where "finite" represents terms that are manifestly well behaved as A„ — 0, and the grouping 
(■ ■ ■) vanishes with A„. Such groupings combine a variety of dilogarithms, logarithms and 
constants together in a non-trivial way. In fact, for higher rank tensor integrals, with higher 
powers of Gram determinants in the denominator, it becomes even more desirable to organise 
the scalar integrals in this way. It is possible to construct these well behaved groupings by 
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brute force, making a Taylor expansion of the scalar integrals in the appropriate limit. 
However, as we will show, they systematically and naturally arise by considering the scalar 
integral in D + 2 or higher dimensions^ and/or by differentiating the scalar integrals with 
respect to the external kinematic variables. Our approach is therefore to re-express the 
formfactor coefficients in terms of functions that are finite as A„ — > 0, explicitly cancelling 
off factors of the determinant where possible. The one-loop matrix elements for physical 
processes will then depend on these finite combinations, which can themselves be expanded 
as a Taylor series to obtain the required numerical precision. An additional improvement is 
that the physical size of the resulting expression is significantly reduced because the scalar 
integrals have been combined to form new, more natural functions. 

Of course, one loop amplitudes may also contain spurious singularities other than those 
directly arising from Gram determinants. Such singularities may occur as one or more of the 
external legs becomes lightlike or as two external momenta become coUinear. Our approach 
has the advantage of avoiding such 'fake' singularities. Since the new finite functions are 
obtained by differentiating the scalar integrals, they cannot contain additional kinematic 
singularities beyond those already present in the scalar integral. This helps to ensure that 
only genuine poles - those allowed at tree level - are explicitly present in the one-loop matrix 
elements. Once again, this helps to reduce the size of the expressions for the amplitudes. 

In this paper, we address how such finite functions might be generated for arbitrary 
processes. We will closely follow the notation and approach of Bern, Kosower and Dixon 
to derive relationships between integrals with polynomials of Feynman parameters in the 
numerators as well as between integrals with fewer parameters but in higher dimensions. 
The basic definitions and notations are introduced in section 2 and the recursive relations 
for integrals with up to four Feynman parameters in the numerator are presented along with 
the dimension shifting relation of ^. These expressions are valid for arbitrary internal 
and external masses and for general kinematics. However, making sense of these relations 
with respect to the singular limit depends on the actual integral itself; i.e. on n and the 
specific values of the kinematic variables. The remainder of our paper describes a series 
of explicit realisations of the three, four and five point integrals relevant for the one-loop 
corrections for the decay of a virtual gauge boson into four massless partons [|1T], |T^. 
In section 3, we consider tensor three point integrals for all internal masses equal to zero, 
but for general external kinematics. Section 4 describes tensor box integrals with one and 
two external massive legs while the tensor pentagon integral is explicitly worked through in 
section 5. Our main results are summarised in section 6, while some explicit results for three 
and four point tensor integrals are collected in the appendices. 



^We note that it will not prove necessary to explicitly compute the scalar integrals in higher dimensions, 
since they will be obtained recursively from the known scalar integrals in Z? = 4 — 2e dimensions. 
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2 General notation and basic results 



The basic integral we wish to work with is the rescaled one-loop integral in D dimensions^, 

r 1 D/2~n 

I^[l] = (-l)"r(n - D/2) / d^x,6il - E S^JX,x, 

•^0 i [i,j=l 

Here, the Feynman parameters Xi have been introduced and the loop momentum has been 
integrated out. The symmetric matrix Sij contains all the process specific kinematics and 
reads, 

. _ {Mf + Mf- (g.^i - 



Following closely the steps of P, O] we perform the projective transformation 



ami 



Xi (XiCli 



En ' 
7 = 1 



i=l 



and introduce the constant matrix pij such that. 



Si' 



Pij 



The parameters ai can be related to the kinematic variables present in the problem, while 
Pij is considered independent of the ai. Provided that all ai are real and positive we find. 



n-D 



Ui 



E Pio^^^i 



D/2-n 



It is useful to rescale the integral. 



(2.1) 



D 

n 1 



SO that in I, the only dependence on the parameters ai lies in the factor J2^=i Differen- 
tiating with respect to a^ brings down a factor of the rescaled Feynman parameter under 
the integral, 

1 di^ii] 



in — D) dai 



(2.2) 



where the notation is obvious. With repeated differentiation, it is possible to generate all 
integrals with Feynman parameters in the numerator. 



^Note that our definition of I differs from that of [pi B[ by a factor of (— 1)*^ 




and, 



The second step of Bern, Kosower and Dixon's work Q is to relate the n-point integral 
with one Feynman parameter in the numerator to a collection of scalar n and {n — l)-point 
integrals, 

In M = ^ E f ^ - ^ J + f-m. (2.3) 

where0, 

ij=l jr=l j=l 

Nn = ^(det 77)^. 

In sections 3, 4 and 5, explicit examples using this notation will be worked through. Equa- 
tion (|2.3|) is the analogue of the formfactor reduction in momentum space of and is easily 
obtained by integration by parts. The summation over m represents all possible pinchings 
of the n-point graph to form {n — l)-point integrals. As expected, we immediately see the 
appearance of the Gram determinants in the denominator. However, equations (|2.2| ) and 
(|2.3| ) are equivalent and since, with a few notable exceptions, the scalar integrals have a 
Taylor expansion around = 0, the act of differentiation will not usually introduce a sin- 
gular behaviour. Therefore, we might expect that the ra-point and (n — l)-point integrals 
combine in such a way that the A„ — > limit is well behaved. We can see how this happens 
by considering the n-point integral in D + 2 dimensions [§, ^ 0] , 



" {n-l-D/2){n-D){n-D-l)^''dajdai 



(n-D-l) A„ \''' ' 2N, 



n m=l 



SO that 



I^N = l)7Jr'[l] - E • (2.5) 

V m=l / 

It is important to note that there are no Gram determinants visible in this equation. They 
have all been collected into the higher dimensional n-point integral. It is clear that if I^[l] 
is finite as A„ — 0, then so is I^"'"^[l] and therefore so is I^[aj]. This confirms that the 
apparent divergence as A„ ^ is fake. Furthermore, I^"'"^[l] is an excellent candidate for a 
finite function - it is well behaved as the Gram determinant vanishes and is easily related to 
the Feynman parameter integrals via equation ( p.4|) . Of course, it may still be divergent as 
e ^ and the dimensionally regulated poles remain to be isolated. 

''Note that our definition of 7 coincides with 7 of ^. 
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By applying the derivative approach, we can easily extend this to two or more Feynman 
parameters in the numerator, 



^ 1 d\^[a. 

'^^ in-D-l) da 



2iV„ Y da, ^ daj " (n - D - 1) ^^"^ 9a, 

Using equation (p^) and the identity. 



we see that. 



= ^ ((^ - ^ - 2)7J^^'K] + - E • (2.6) 

Note that ifi™''[l] does not depend on am, and therefore, 

p(j)r ^ _ ^In-y[l] _ ^ 
In-1 K] ~ = 0. 

Consequently, the m = j term in the summation vanishes. 

Differentiation has not produced any new Gram determinants and we can treat these 
integrals as new well behaved building blocks, or substitute for them using equation (p.5|) 
with D replaced by D + 2, 

InW = ^((n-/^-2)(n-D-3)Wf+1l] + 2Ar„r^,,if^ 



-{n-D-2)Y: l^V,Jn-l^"'\l] - 2iV„ Y: V^Jn^T\c 
m=l m=l 



The scalar integrals for D + A dimensions can be obtained recursively from equation (|2.4|). 



Replacing the factors of a^ in equations (|2.5|) and (|2.(j|) and the analogous equations for 



three and four Feynman parameters in the numerator, we find, 

InN = ^((^-^-l)7Jf+'[l] -E^.^«,amlfir^[l]), (2.7) 



In k.^.] = ^ ((^ - ^ - 2)7Jr'[a;,] + r/.,a,a,lf +2[1] 



- Vimaiamln-T\^j] ) > (2-8) 

m=l 



[xiXjXk] = \{n- D- 3)7il^+^ [xjXk] + rjijaiajl^^'^ [xt] + r/ifcttiafcl^^^ [x^] 



n 

fD(m) 

m=l 



^ ^ VimCXi(^mJ"n—l [-^j-^k] j ; (2-9) 



[xiXj-ajfcXi] = Un- D - 4)7J^+^ [xjXfcX/] + r]ijaiajl^^'^ [xkXi] + r]ikaiakln^^ [xjXi] 



+r]iiaiail^^^[xjXk] - Vimaiamln-T\xjXkXi] ]. (2.10) 

m=l 



Once again, no Gram determinants are apparent and these equations may be solved by 
recursive iteration. These are our main results and their use will be made clear with the 
explicit examples in the following sections. 

Before proceeding to the explicit examples, we note that the full tensor structure in 
momentum space is simply obtained from the Feynman parameter integrals by introducing 
the auxiliary momentum defined in equation ( |1.3| ). With an obvious notation (and after 
integration of the loop momentum) the tensor integrals can be written, 

where {a. . , b^t^^-t^" is the usual Passarino-Veltman notation 0, and indicates a sum over 
all possible permutations of Lorentz indices carried by a ... 6. For example. 



Throughout the next sections, we make the simplifying choice that Mj = 0. Such integrals 
are relevant for a wide range of QCD processes involving loops of gluons or massless quarks. 
The approach can be straightforwardly extended to include non-zero masses P]. 

The strategy is to isolate the ultraviolet and infrared poles from the tensor integrals, 
leaving the finite remainder in the form of groups of terms that are well behaved in all of 
the kinematic limits. In real calculations where groups of tensor integrals are combined, this 
grouping will often cancel as a whole. Alternatively, if the kinematic coefficient allows, the 
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determinant can be cancelled off for all of the terms in the function. This approach is well 
suited to treatment by an algebraic manipulation program, once the raw integrals have been 
massaged to isolate the poles in e and to group the terms. As we will show in the explicit 
examples, this is usually straightforward. 



3 Three point integrals 



In processes where the internal lines are massless, there are only three types of triangle graph 
described by the number of massive external legs. For the one loop corrections to five parton 



scattering ||T5|, |T6|, ^Tj, only the graphs with one and two massless legs occur. For processes 
involving a gauge boson such as Z — > 4 partons []1T], [1^, |13|, graphs with all external legs 
massive or off-shell contribute. 



3.1 The three-mass triangle 

We first consider triangle integrals with exiting momenta pi, p2 and ps as shown in fig. ^ 
and all internal masses equal to zero, Mj = 0. Throughout, we systematically eliminate ps 
(and the Feynman parameter X2) using momentum conservation so that p^ = —{pi +P2), 
pI = {Pi+P2f = S12 and, 

V>' = -{l-x,)p't~Xspi^. 

The full tensor structure with up to three loop momenta in the numerator can therefore be 
derived from loop integrals with up to three powers of xi or X3 in the numerator. 

As a first step, we consider the general case, pf, pi, su 7^ 0, where the scalar integral in 
four dimensions is known to be finite. Here the ctj parameters can be determined by. 



while 



A3 = -Pt -Pt- Su + '^pIpI + 2j5iSi2 + 2p^Si2, 

A3 = —a\ — a\ — a\ + 2a\(X2 + 2a\0iz + 2020:3. 

From the definition of the matrix 77, we see that, 

/ -1 1 
= 1 -1 1 I , A^3 
111 
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«2 



P1 +P2 




«2 



Figure 1: The triangle graph and each of the three pinchings obtained by omitting the 
internal line associated with am for m = 1, 2 and 3. In each case, the internal line is shrunk 
to a point and the momenta at either end are combined. The relation between the external 
momenta and the a, can be seen by cutting the loop; aittj = — where p is the momentum 
on one side of the cut and at, aj are associated with the cut lines. 



The variables 7 always appear in the following combinations, 

^. ( pUpI + S12 - pI) « = i 
' 2 = { ^i2bi + pI- S12) i = 2 

n.=i«. [ p2(^p2 + - pI) ^ = 3 



The scalar triangle integral for all external masses non-zero is finite in four dimensions 

m i i n 



= (^log(a+a-) log ([3^) + 2Li2(a+) - 2Li2(a-) j , (3.1^ 

where Li2 is the usual dilogarithm function and are two roots of a quadratic equation. 



a 



± +pI-pI± 



2Sl2 



Although l3[l] appears to diverge as A3 — > 0, this is not the case. As noted by Stuart |T0 
in this limit, the triangle graph reduces to a sum of bubble graphs, 

lim Isfl] = % ^log H % ^log ( ^ I , 

and there is a well behaved Taylor series in A3. 

^For scalar integrals in D = 4 — 2e or = 4 dimensions, we omit the superscript D. 
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3.1.1 Tensor integrals in D = 4 



The tensor integral can be easily written in terms of higher dimensional scalar integrals and 
bubble scalar integrals using eqs. (|2.7| - p.9|) . For one Feynman parameter in the numerator, 
this gives, 

^ m=l 

Immediately a problem is apparent - the coefficient of the scalar integral in (6 — 2e) dimen- 
sions, 7j, is singular as one or more of the external momenta become lightlike. Although the 
divergence as the Gram determinant vanishes has been removed, it appears to have been re- 
placed by a divergence as the invariants vanish^. However, these divergences cancel between 
the triangle and bubble contributions and the tensor integral itself is well behaved and finite 
in all kinematic limits and is therefore a better choice for a finite function. 

In fact, since the triangle scalar integral is finite in 4 dimensions, it is convenient to 
generate the tensor structure directly from derivatives of the scalar integral. However, in 
order to use equation ( p.3| ), we also need the two point integral for external momentum p 
(and internal masses Mj = 0), 

for each of the three pinchings, m = 1, 2 and 3 shown in fig. For pinching m and D = 4—2e, 



1-e 



^UUc^^J e(l-2e) Vnti«i/ 

where the usual product of Gamma functions obtained in one-loop integrals cr is given by, 

r^(i-6)r(i + 6) 
^ r(i-2e) ■ 

Rewriting equation ( |2.3| ) for the case D = 4, n = 3 and i = 3 and adding, 

2^=iV A3 / a2 



we see, 



Has] = i f73i3[l] + I E (737.. - r^3^A3) ^ (ir^[l] - if [1] 

^3 \ ^ m=l 11" 



^ '^3i3[l]-2iiog(^^)+21og(^))+0(6). (3.3) 



Aq V "3 Vai/ Vas 



'Problems in this limit are to be expected since even the scalar integral itself is not finite as pf 0. 
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Alternatively, this could be obtained by differentiation of equation (|3.1|) . By trivial replace- 
ment of factors of a, we find, 



^ U(^i2 + pI- P\)h[^] + {P\ +pI- ^12) log te j - 2p\ log f ^ 



Integrals with higher powers of Feynman parameters can now be generated by direct differ- 
entiation of I3[a3], 

1 d 

h[x^X^^] = -(^^^ (n«) «^«3^l3[a3]- 

All of these functions will be finite in the e ^ limit, and can be considered as building 
blocks in constructing the tensor structures for box and pentagon integrals. In fact, because 
they are obtained by differentiating a function well behaved as A3 0, they are also finite 
in this limit. Therefore, they tie together the dilogarithms from the triangle integrals and 
the logarithms from bubble integrals in an economical and numerically very stable way. 

We also see that they are directly generated in tensor structures for box graphs (equa- 
tions ( p.7| - |2.1(i| ) with n = 4) and will naturally cancel in Feynman diagram calculations 
involving both triangle and box graphs. For general calculations with pf ^ and pi 7^ 0, we 
introduce the notation, 

Lco(pi,P2) = Isll], Lc2„_i(pi,P2) = hl^s], Lc2n(pi,P2) = Mxi^s], (3.5) 

for LC0...5. The symmetry properties of the triangle function imply that the analogous func- 
tions for xi ^ X3 (or ai ^ as) are just obtained by exchanging pi and p2- In dealing 
with box graphs, integrals with X2 in the numerator will naturally arise. In these cases, we 
systematically eliminate them using J^i^i = Explicitly, we find, 

hlxixs] = ^^(^plisu + pI - pl)h[x3] + pI{si2 + pI - pl)h[xi] 

-py,h[l]-pllog(^j +P?+P2-Si2j, (3.6) 



^3] = ^f3p?(.i2+P^-p?)l3N+Pll3[l]-(^12-p^)l0gtej-2p? 



I3[a;ia;3] = y^pl{s 12 + pj - pl)h[xl] + Qpl{s 12 + pI - pl)h[x 1X3] 

-Splplhlxs] + Pih[xi] - pI log +pI- si2^ , (3.8) 

h[xl] = ^f5p?(.i2+P^-pDl3[a:^] + 2ptl3[x3]-(3i2-p^)logfe') -p?V(3.9) 
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10''^ 10" 10-" 10^'' 10"" 10"' 10''' 10'^ 10"* 10"^ 10'^ 10"' 10° 



Figure 2: The finite functions for the triply massive triangle graph with Si2 = 1, Pi = 0.2 as 
a function of As/A™^^ where A™^^ = — (si2 —pf)"^. The functions have been evaluated using 
double precision Fortran. The dashed lines show the approximate form for the function in 
the limit A3 — > 0, retaining only the first term of the Taylor expansion. 

By expanding as a series in A3, these functions can be evaluated near the singularity 
with arbitrary precision. For example, 



To illustrate this, fig. H shows the various functions at a specific phase space point, 
S12 = 1, pI = 0.2 and letting vary in such a way that A3 — 0. This corresponds 
to P2 0.135. We see that as this limit is approached, the numerical evaluation of the func- 
tion using double precision (an intrinsic numerical precision, acc, of roughly 10"^"^) becomes 
uncertain. For this particular phase space point, functions with a single Gram determinant in 
the denominator (Lci) remain stable until A3/A™'^^ ~ 10~^ while those with more powers of 
the Gram determinant break down correspondingly sooner - at A^/ A'^^^ ~ 10~^ for Lc2 and 
Lc3 and A3/A™'^'^ ~ 10~^ for LC4 and LC5. In general, numerical problems typically occur 
when A3/A™^^ ~ (acc)^/^ where N is the number of Gram determinants in the denominator 
of the function. Other phase space points yield a similar behaviour. 

The unstable points represent a rather small proportion of the allowed phase space. 
However, problems may arise using adaptive Monte Carlo methods such as VEGAS 






3si2(si2 -pj- pI) 3si2(si2 - Pi - plf 



(3gi2-3pf ^(si2\ 

6si2(si2 -p? -pi) KpI)' 
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where the phase space is preferentially sampled where the matrix elements are large. Finding 
an anomalously high value for the matrix elements in a region of instability would cause the 
Monte Carlo integration to focus on that region giving unpredictable results. 

Of course, these instabilities could be handled by a brute force increase in numerical 
precision. While possible, this has the disadvantage of producing significantly slower code, 
and, since in all cases, the approximate form obtained by making a Taylor expansion about 
A3 = and keeping only the constant term works well where the numerical instabilities begin, 
this is not an attractive solution. In fact, the approximation is reliable for A3 < 10~^A™'^^. 
Explicit forms for the approximations are collected in Appendix C. 



3.1.2 Scalar integrals in higher dimensions 

We now turn to the scalar triangle integrals in higher dimensions. They appear in the g^v 
part of the general Lorentz structure and recursively in the determination of if* 
Unlike the triangle in four dimensions, these integrals are ultraviolet divergent due to the 
presence of the various pinchings - bubble integrals. A function that can usefully be used as 
a building block of matrix element calculations, must be finite as both A3 — >■ and e — >■ 
and we must first isolate the poles in e. Although equation ( |2.4| ) suggests that the ultraviolet 
pole structure involves A3, this is easily shown not to be the case. Adding terms proportional 
to, 



1 f A3 - E 7n^«-l -f^[l]=0. 



to equation (|2.4| ) for n = 3, D = 4 — 2e we see. 



where the divergence as e ^ lies exclusively in the last term. Reinserting the factors of a 
and the definition of l2[l] in D = 4 — 2e dimensions we find. 



l3^=^-^^[l] = Lci,(pi,p2) - \ + 3J cr, (3.10) 

where, 

-^2(^12 + Pi - vl) log ) 



— e 



Up\h[xi]+plh[x^]). (3.11) 
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In a similar fashion, the e pole structure can be removed from the triangle scalar integral 
in D = 8 — 2e and D = 10 — 2e dimensions yielding two more functions that are finite as 
both A3 ^ and e ^ 0. Explicitly we find, 

r-ll] = W(p..p,)-M±|±£l2)((z£l2n + |j,„ (3.12) 

(3.13) 

where the finite functions are defined by, 

Lc25(pi,P2) = ^ (^2pIpIsi2Lcis{puP2) - ^ (^pt{s 12 + pI - pI) log 

+Pt{si2+pl~pl) log f ^ ) + 2pIpIsu]] , (3.14) 



P2 

LC35(Pl,P2) = ^^hplplsi2^C2s{Pl,P2) - ^(pt{Sl2+pl-pl)^Ogl^ 



+Pt{si2 +pI- pI) log + +pI + s,2)\ \ . (3.15) 



Although these functions have been obtained via equation (|2.4|) , they are still related to 
derivatives of the basic scalar integral in D = 4 — 2e, and are therefore finite in the A3 ^ 
limit. We can see this by examining the same phase space point as before, S12 = 1, pi = 0.2 
and varying A3. As expected, the Lc„5 show a similar behaviour to the Lc,i functions 
- numerically breaking down at larger and larger values of A3 as the number of Gram 
determinants increases, and being well described by the first term in the Taylor expansion 
as this happens. For completeness, we collect the limiting approximations in Appendix C. 



3.1.3 Tensor integrals in higher dimensions 

For triangle loop integrals with three loop momenta in the numerator, it is also necessary to 
know the D = Q — 2e integral with a single Feynman parameter in the numerator. Rather 
than differentiating the ultraviolet divergent lf~^~^'^[l], we can evaluate it in terms of the 
D = A — 2e tensor integrals of section 3.1.1. Using ( |2.6| ) for D = 4 — 2e, we see that, 

h[a,a,] + i3[a3a,] = \ (^-(3 - 2e)(7i + %)l^='~^'[a,] + (r/i, + 

- Z!(^l™ + ^3m)l2'"^K]J, 
m=l / 
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which, for j = 1, 3,[] simphfies using, 

vij + mj = 0. 

The same equation simphfies the sum over bubble pinchings so that only m = 2 contributes, 
while 7i + 73 = 2a2- Restoring the factors of a and using. 



yields. 



-2er^ 1 _ 1 / 2t r 1 , ^2t r n \ 1 / (-■5l2, 



— e 



ir'l^i] = ^[pih[xix,] +pth[x3X,] - - [ I + 3 cr. (3.16) 



Later, we will see that constructing the tensor integrals for box graphs can also generate 
l3~^^[xj] and higher dimensional triangle integrals with more parameters in the numerator. 
In each case, we use similar tricks with equations ( |2.7| - |2.lGD to rewrite them in terms of the 



four dimensional integrals with the ultraviolet pole made explicit. 



3.2 The two-mass triangle 

We will also be interested in triangle graphs where one or more of the external momenta is 
lightlike. Here, we first focus on the — * 0. In D = A — 2e dimensions, we have the 
well known result, 

r[i] = ^f^^i^^ii^i^V (3.17) 

V "3 - "2 / 

where the superscript indicates that only two of the three legs are massive. For this choice 
of kinematics, 

aia2pl = —1, aia3Si2 = —1, 

and 

A3 = -(03 - a2)^ 

so that the singular limit is 03 — >■ 02. Because A3 makes no reference to ai, rjij contains a 
row of zeroes, 


-1 
1 



Vij 



and therefore N3 = 0. Consequently, care is needed in applying the equations of section 2. 
In addition, since the scalar integral for three massive legs is finite (and the results in the 
preceding subsections have been explicitly derived in D = 4), one cannot just set P2 0. 

^Since the Feynman parameters add to one, the case j = 2 is of httle interest. 
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In fact, it is easiest just to bypass the problem and generate the whole tensor structure 
by direct differentiation of the scalar integral with respect to 0:3 and ai. It is easy to see 
that, 

I^K] ~ + B^^^ + Clog (^) + D, 

where A,B,C and D are polynomials in 1/0:3 and 1/(0:3 — 02). Unlike the all massive 
case discussed before, the scalar integral is singular as e — > 0. As a general rule it is not 
necessary to be particularly careful with double poles in e, since they must either cancel 
or form the infrared poles of real matrix elements. However, it is possible for the integrals 
to be multiplied by factors of e - from expanding factors of dimension - and the resulting 
logarithms should occur in combinations that are finite as 03 ^ 02- It is easy to see that, 

log f ^) 

03 — 02 

is finite. So, to tie the logarithms and constants together in combinations that are well 
behaved in the 03 — 02 limit, we use the fact that derivatives of this function are also well 
behaved, and introduce the functions, 

Lcr(pi,P2) = - Jim (e X iffxr^]) , (3.18) 
for n = 1 4. In terms of invariants, 



T 2m/ N /'PlLc^-l(Pl'P2) — ;^ , /o in\ 

Lcr(Pi,P2) = - I 2 — I ' ^ ^ 2 (3.19) 

■S12 — Pi 



with. 



Lc^(pi,P2) = Lcr(Pi,P2) = log ('-^-] . (3.20) 

S12 -Pi \Pi J 

These functions, or functions closely related to them, have appeared in next-to-leading order 
matrix element calculations [|1^, |1^, |1^, The explicit forms for l|'"[a;iXj] appearing in the 
momentum expansion are well known and are collected in Appendix A. 

Although these functions are rather simple, they still contain numerical instabilities as 
Pi Si2. This can be seen in Fig. | where we show Lc^™" for the specific phase space point 
S12 = 1 and let pf approach S12. While a single inverse powers of {s 12— pi) is handled correctly, 
higher powers cause problems. As can be seen from the figure, a suitable approximation is 
obtained by the first term in the Taylor expansion, 

Jim Lcl"'{pi,p2] 



pf^si2 npi 

for n > 1. 
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10'' 10* 10"' 10'" 
(Si2-pf)/Si2 



Figure 3: The finite functions for the triangle graph with two external masses with Si2 = 1 
evaluated in double precision Fortran as a function of (si2 —pD/su- The dashed lines show 
the approximate form for the function in the limit pi — > Si2, retaining only the first term of 
the Taylor expansion. 



The other configuration of triangle graph that appears is where two of the momenta are 
lightlike, Pi = P2 = 0. Once again, the tensor structure can be generated by differentiation 
or canonical Passarino-Veltman reduction. Here, there is only one scale in the problem 
so there can be no logarithms and it is neither possible nor necessary to introduce well 
behaved functions. The explicit forms for the Feynman parameter integrals appearing in the 
momentum expansion are well known and for the sake of completeness are given in Appendix 
A. 

This concludes our discussion of triangle graphs. For the case of three massive external 
legs (and internal masses set equal to zero) the four dimensional tensor integrals are finite 



as e — > and are given by the functions Lcq... 5(^1,^2) defined in ( |3.1| , p^p^6| - p.9|) , while the 
ultraviolet divergent part 5^^,^ part, if "^"^^[1] is expressed in terms of a similar function (Lci^) 



with the pole isolated ( |3.11|) . For the tensor structure in the simpler case with one lightlike 
leg, it is useful to group the logarithms and constants using the functions Lcl"^4^{pi,p2) ( ^.2(J| , 
TT9D. 



For the more general case where the internal masses are non-zero, the same procedure can 
be utilised. The matrix rj has slightly more entries and there are more scales in the problem. 
However, the grouping together of triangle graphs and bubble integrals into functions well 
behaved in the A3 — limit and the isolation of the ultraviolet singularities can be made 
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Figure 4: The box graph and each of the four pinchings obtained by omitting the internal 
hne associated with for m = 1, 2, 3 and 4. 

exphcit in the same way. 



4 Four point integrals 

For one-loop corrections to five parton scattering, box graphs with at most one external leg 
occur. However, for processes involving a massive vector boson and four massless partons, 
we can obtain box graphs with a second massive external leg by pinching together two of 
the partons. There are two distinct configurations according to the positions of the massive 
legs; the adjacent box graph and the opposite box graph. The box graph is shown in 
fig. ^ for outgoing momenta pi, p2 and p^. Throughout this section, we will assume that 
(Pi + P2 + PsY = ■Si23 7^ 0. In the adjacent two mass case, pi = pi = and pf ^ 0, while 
for the opposite box, Pi = pi = and p^ 0. Unfortunately, the raw scalar integrals 
for these two cases behave rather differently. The adjacent box is finite in the limit that 
A4 — >• 0, while the opposite box diverges as A4 — > 0. In this section, we work through these 
two configurations and rewrite the tensor integrals in terms of well behaved functions and 
explicit poles in e. 
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4.1 The adjacent two-mass box 



We first consider the adjacent box with P2 = P3 = and all internal masses equal to zero. As 



in the triangle case, we systematically eliminate one of the momenta, 
and one of the Feynman parameters, X3 = 1 — xi — X2 — X4, so that. 



'p^' 



[1 - Xi)pl2 + X2P^ - Xip'^. 



The related integrals with pi ^ and pf = are obtained by pi ^ p^ (and the indices i and 
j in and rjij transform as 1 4 and 2^3). 

For this kinematic configuration, the a parameters are defined by, 

^2 



-1, aia2pi = -1 
-1, 0204523 = -!, 



while, 



Vij 



( 


1 

V 






-1 
1 



1 0\ 

-1 1 

2 -1 
-1 ) 



Na 



1 

2' 



and, 



A4 = 2S23 ((Sl23 - Si2)(Si2 - p\) - S12S23) 



The coefficients 7^ always appear in the following combinations which are directly related to 
the conventional variables, 



It 



n 



3=1 



-S23 

S123 - Si2 
S123 +pI- S23 
I pI - Sl2 



S12 



1 

2 
3 
4 



The scalar integral in D = 4 — 2e can be written 



cr 



14(1] = — ((a2a4)' + 2(aia3)' - (aia2)' - {aiaAY) + 2Ldo(pi,P2,P3) + 0{e), 



(4.1) 



where. 



Ld„(p„p„P3) = Li. (1 - I) - Li, (1 - I) + ilog (i^) log (I) . (4.2) 



In constructing the tensor integrals in D = 4 — 2e, we see from equations ( |2.7| - |2.10D that 
the box integral in higher dimensions is needed. In fact, in D = 6 — 2e, the box integral is 
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infrared and ultraviolet finite. This can be seen by inspection of equation (|2.4| ) and noting 
that the pinchings with m = 1, 2 and 4 in the expression, 

m=l 

are proportional to and, when combined with the appropriate 7 factor, precisely cancel 
with the pole structure of the box integral. The final pinching (m = 3) corresponds to the 
triangle graph with three massive external legs which is itself finite. Altogether, we find that 
the adjacent box integral in D = 6 is, 

I4 [IJ = I Ldo(pi,P2,P3) + 2 hl23 +Pl - S23 I Lco(pi,P23) I 

= Ldi5(pi,p2,P3), (4.3) 

where Lco(pi,P23) = l3[l] is defined in equation Because of the finiteness properties of 

the three mass triangle, we will find repeatedly that the m = 3 pinching should be treated 
differently from the other three. 



4.1.1 Scalar integrals in higher dimensions 



For higher dimensions, we just reuse equation (|2.4|) , noting that the triangle pinchings in 
D = 6 — 2e reintroduce ultraviolet poles. These can easily be isolated by adding and 
subtracting combinations of scalar integrals as in section 3.1.2. Explicitly, 

= L4,fe,p„p3) + |(t£2l)2 + ^j, (4.4) 

= Ld3.(p..p..P3) + + + + ^) ■ (4,5) 

HI ^ Ld„(p.p.p3) + |^(<^.f). (4.6) 



where, 



9 9 9 A 9 9 9 ^12^23 

-P = ■S123 + ■5l2 + -523 + Pi + -3123512 + S123S23 + S123P1 + ■5l2Pi + S23P1 H 7, ■ 



The finite parts of the higher dimension boxes are given by, 

Ld25(Pl,P2,P3) = -^^^■Sl2S23Ldi5(pi,P2,P3) 
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+ Si23 +Pi - S23 Lci5(pi,p23j 

2 / 

523, /Si23\ , /■5i23\ Pi, / ^123 




Ld35(pi,P2,P3) = -^^^|^Si2S23Ld25(Pl,P2,P3) 

+ S123 +Pi - S23 LC25(Pl,P23 

V ^12 / 

-,2 /„ N „2 



+^ log f ^) + ^ log V - log + ^ I (4-8) 

24 ^ V S23 / 12 ^ V S12 / 24 ^ I pf y 12 r ^ ^ 



Ld4s(Pl,P2,P3) = -^l^^f Sl2S23Ld35(pi,p2,P3) 



+ S123 +Pi - S23 Lc35(pi,p23) 

V ^12 / 

log f^U 4 log - t4 log { '-^ 



360 V S23 / 180 V S12 / 360 \pi 

, Si2S23(si23 + S12 + S23 + pI) \ 

+ 720 )■ ^^-^^ 

The D = 6 — 2e and D = % — le box integrals explicitly appear in the momentum space 
tensor structure with one and two factors of g^y respectively. All of these integrals appear 
either directly or indirectly in the tensor box integrals of equations (|2.71 - |2.ld|) . 

Once again, all of these functions are well behaved as A4 and group a variety 
of dilogarithms, logarithms and constants together in a non-trivial way. This is shown in 
Fig. H for a particular point in phase space; S123 = 1, S12 = 0.4, S23 = 0.08 with p{ varying 
so the A4 ^ limit is approached. We see that although Ldis', with a single inverse 
power of the Gram determinant, is numerically stable, the functions with more powers of 
Gram determinant in the denominator break down at much larger values of A4. In all 
cases, the function is well approximated by the first term of the Taylor expansion provided 
A4 < 10~^A™'^^. These approximations are collected in Appendix C. 



4.1.2 Tensor integrals 

Armed with the integrals in D > 4, we return to the tensor integrals and use equations ( p.7| - 
p.lOl) as the starting point. Rewriting equation (|2.7|) for n = 4, D = 4 — 2e and noting that 
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10 '^ 10 " lO-'" 10-^ 10-' 10-' 10-^ 10-= lO* 10-^ 10-2 10-' 10" 



Figure 5: The finite functions for the one mass box graph as a function of A4/A™'^ where 
= 2si2S23(si23 - Si2 - 823)- The phase space point is S123 = 1, 512 = 0.4, S23 = 0.08 

and p1 altered so the hmit is approached and the functions have been evaluated in double 
precision Fortran. The dashed lines show the approximate form for the function in the limit 
A4 — > 0, retaining only the first term of the Taylor expansion as given in Appendix C. 
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we have eliminated so that i = 1, 2 and 4 only, we have, 

4 



m=l 



The factor ajOm multiplying the triangle pinchings will always produce a factor of For 
triangle graphs with at least one massless leg (pinchings m = 1,2 and 4), the contribution 
is ~ 1/e^ and will combine with similar poles from other Feynman diagrams. On the other 
hand, the m = 3 pinching (corresponding to the triangle graph with momenta pi and P23 
flowing outwards) is finite and, provided the value of s is related kinematically to that 
triangle pinching, there may be a possibility of cancellation with other triangle Feynman 
graphs. However, for the case i = 1 and m = 3, the associated invariant mass is S12. 
This term cannot combine with any other naturally generated triangle graph with the same 
kinematics. Therefore, we group this term with the box integral, by adding and subtracting. 



m=l 



SO that, 



with. 



h[xi] = 7iLdi(pi,p2,P3) - J2 Vimaiam {^3^^^ - ^^i^) , (4-10) 



m=l 



Ldi(pi,p2,P3) = - (lf=ni] + lf [1]) . (4.11 



The only non-zero entries in rjim for i ^ 3 and m 7^ 3 are i = 2, m = 4orz = 4, m = 2 
corresponding to aiOm = — I/S23 which is appropriate for ^'^^[l]. 

For integrals with more Feynman parameters it is convenient to introduce the following 
functions, 

Ld„ii„,i„_i(pi,p2,P3) = -\nlf='^[xi,...Xi^_,] + l'i\xi,...Xi„_,]\, (4.12) 



for n = 2,3 and 4. Suppressing the arguments of Ld^, and using equations (p.8[ - |2.10| ), we 
find, 

h[xiXj] = 7iLd2j - VijO^iOij] Ldi + [1] f - XI VimOiiaml Is™^ [xj] - if ^ [xj] I, (4.13) 
h[xiXjXk] = jiLdsjk - ^r]ijaiajl^Ld2k + ^Pixklj - ^r]ikaiak^Ld2j + l3\^j]^ 

- Y Vimaiam \ Is""^ [xjXk] - [xjXk] > , (4. 14) 

m=l ^ 
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-^Viio^iaiS^Ldsjk + lf\xjXk]^ - ^ r]i^aiam^l3^\xjXkXi] - I'i^XjXkXi]^. 

(4.15) 

Since we have systematically eliminated X3 using the delta function, i, j, k and / run 
over 1, 2 and 4. This guarantees that the coefficients of the form rjijaiaj are only non- 
zero for i = 2 and j = 4 (or vice-versa). In these cases, ?724a2tt4 = —^/s23, which is again 
appropriate for the m = 3 pinching to form a completely massive triangle, ■* [1] . Altogether, 
equations ( |4.10| , [4.13H4.15| ) are sufficient to completely describe the tensor structure of the 



adjacent box. 

However, in order to determine the Ld„j^,..j„_-^ combinations, we need tensor integrals 
for D = 6 — 2e dimensional box graphs with two or more Feynman parameters. These can 
be obtained from equation ( [2.8[ ) once the D = 6 — 2e box integral with a single Feynman 
parameter, lf^^~^'^[xj] is known. This can be derived by differentiating ( p.2|) which indicates 
that is also finite as e — 0. To see this, we reuse equation (p.7|) and our usual 

trick of adding and subtracting combinations of the m = 3 pinching in D = 6 — 2e, 



lf='-''[^^] = 7.(-(3-26)lf=«-2ni]-ir'~'^^'^[l] 



m=l 



Both brackets are separately finite. First, the divergent part of if ^ ^^[1] precisely cancels 
against that of I^~^~^'^*^^'' [1]. Second, all triangles in D = 6 — 2e dimensions have 1/e poles 
and the difference of any two, is either zero or a log. Further differentiation does not change 
the finiteness properties of the lf^^~^^ tensor integrals. 

The Ld„i^...j^_j combinations are also well behaved in certain kinematic limits. For ex- 
ample, If='^[l] and if [1] are finite as S12 or S23 — * 0. Just as differentiating functions 
which are finite as e ^ does not introduce poles in e, neither can it introduce poles in the 
kinematic invariants S12 or S23. As an example, consider the function Ld24 given by, 

Ld24(]?l,]?2,]33) = — (3Ld25(Pl,]32,P3) + LCis(pi, P23)) - — Lci(p3,Pi2) 

S12S23 S23 

+ — Lci(p23,Pl), 

■S23 

which appears to contain a pole in Si2. In the S12 — > limit, A4 — > — 2s23Si23j9i and 
Ldas ^^^|^^Lcis(pi,P23) -^Lci5(pi,p23), 
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so that, 

lim (3Ld25(pi,P2,P3) + Lci5(pi, P23)) ^ 0, 

and therefore, 

hm S12 X Ld24(pi,P2,P3) 0. 

Similarly, Ld24 contains no power-like divergence^ in the S23 — > limit and, with a little 
more work, it can be shown that, 

lim S23 X Ld24(j9l,J92,P3) 0. 

S23^0 

Once again, these functions combine dilogarithms, logarithms and constants in a highly non- 
trivial way to form well behaved building blocks. Explicit forms for the Ldni^...i^_^ functions 
for n = 1, 2 and 3 are given in Appendix B. 



4.2 The one-mass box 



The higher dimension and Feynman parameter integrals for the one-mass box integral ob- 
tained by taking pf can also be constructed in a similar way. For this kinematic 
configuration, the a parameters are defined by. 



aia4 S123 



-1, 



aias S12 



-1, 



02^4 -523 



-1, 



while. 



and, 



In terms of invariants. 



Vij 



( 


1 
V 






-1 
1 



1 0^ 
-1 1 

00/ 



1 

2' 



, Im 



n 



2S12S23 (Si23 — S12 — S23) . 



-S23 i = 1 

S123 - S12 « = 2 

S123 - -523 ^ = 3 

-S12 i = 4 



^ Although these functions do not behave as inverse powers of the vanishing kinematic variables, they do 
contain logarithms of S12 and S23. This is because the e ^ limit has already been taken, and the order of 
taking the two limits does not commute. For next-to-leading order calculations, we only approach the limit 
Sij —^ and e can safely be taken to zero first. 



26 



The scalar integral can be written, 



= ^ {2{aia;r + 2(«2«4)^ - 2{a^a,r) + 2Ldr(pi,P2,P3) + 0{e), (4.16) 



where, 

Ldr(,,p.p3)^L.(l-|).L,(l-^)..og(^).„,(|)-^, 

As expected, in Z) = 6 — 2e dimensions, the scalar integral is finite. 



I^^=^[l] = -%5^Ldr(Pi,P2,P3) = Ld}-(pi,p2,P3). (4.18) 



In higher dimensions, the scalar integrals satisfy analagous equations to ( |4.4| - |4.6| ) with 
— > 0, and finite parts given by, 



Ld25(pi,P2,P3) = -|^l^^Sl2S23Ld}5(pi,p2,P3)+S23log(^^^ +Si2log(^^^^^ 



(4.19) 



Ld^5(Pl,P2,P3) = -^(^12^23Ld^5(Pl,P2,P3) + ^l0g 



/Sl23\ 


s2 

+ — log 




V S23 / 


12 ^ 


V Si2 / 



+^), (4.20) 



^4"' V i»U \ S23 / i»U \ Si2 

i3 + 

720 



^£H£H(£2±£ll±iH) . (4.21) 



The stability of these functions as A^™ — > is illustrated in fig. ^ for a particular point in 
phase space; S123 = 1, S12 = 0.3, with S23 varying so the limit is approached. The maximum 
possible value of A^™ occurs when S23 = (S123 — Si2)/2; i.e. A4" ^^"^ = Si2(si23 — Si2)^/2. As 
before, Ldj^™ , with a single inverse power of the Gram determinant, is numerically stable. 
However there are numerical instabilities for the other functions with more powers of Gram 
determinant in the denominator. In all cases, the function is well approximated by the first 
term of the Taylor expansion provided A4™' < 10~'^A4™ 

Unlike the pf ^ case, the scalar triangle pinchings all contain infrared poles and there 
is no benefit in absorbing the m = 3 piece in the tensor integrals. Therefore we introduce, 

Ld^™...i„_i(Pi>P2,P3) = -nl4'"'^='^[xi, . . .Xi„_J, (4.22) 
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Figure 6: The finite functions for the box graph with two adjacent massive legs as a function 
of A\"^/A\"^ The functions have been evaluated in double precision Fortran. The dashed 
lines show the approximate form for the function in the limit A\"^ 0, retaining only the 
first term of the Taylor expansion as given in Appendix C. 



for n = 1, 2, 3 and 4. Using equations ( |2.7| - f^.lOD with n = 4 and D = 4 — 2e we find 



1-4 



m=l 



m=l 
1 



Im 



A 

im) r 



E-i\m) r 



m=l 



(4.23) 
(4.24) 



(4.25) 



111 



r]imOiiaml'3'\xjXkXl]. 



(4.26) 



m=l 



As in the previous section, the Ld;^™ functions are finite as e ^ and contain no 
power-like divergences in the Si2 — > 0, S23 — > and A\"^ — > limits. For convenience, explicit 
forms are given in Appendix B. 
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4.3 The opposite two-mass box 



The two mass box graph where the massive legs sit on opposite sides is a special case because 
the scalar integral itself is not finite as 0. We must therefore proceed with care. To 

make best use of the symmetry under pi pa, it is convenient to write, 

= -(1 - xi)p'^ - (x3 + a:4)p2 - 2^4P3- 
Under this flip symmetry, Xi x^, X2 x^ and 



In this kinematic configuration, the a parameters can be defined by, 

ttia4'Si23 = -A, a2a^p'i = —l 
ttiQis S12 = —1, 0.20.4^ S23 = —1, 

where A is an extra kinematic variable that ensures that the are independent, 

A 



Sl23Pi 



With this choice of a,. 



A) 



/ 








1 


-1 \ 










-A 


1 




1 


-A 








\ 


-1 


1 





) 



Each row of rj naturally couples together two of the pinchings (triangles) of this box, so we 
might expect such structure to dominate the integrals. The associated Gram determinant is 
given by. 



^opp _ 2(5^2523 ~ P2'^123)'5l3 



2(1-A)si3 



n 



i=l 



where, 
and. 



Si3 = S123 - S12 - S23+P2, 



A) 



P2 - S23 ^ 

•5l23 ~ 'S12 i 

S123 - S23 i 



1 

2 
3 

4 



P2 - S12 I 

We note that the presence of the Gram determinant is synonymous with a factor of (1 — A). 
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The scalar integral for the opposite box in D = 4 — 2e is given by, 



ir[l] = (y4a) (5 (^"'"^^^^ + ~ ~ («i«4)^A-^) + Ldr(Pi,P2,P3))+0(e), 

(4.27) 

where the finite part Ldo can be written, 

Ldr(pi,P2,P3) = Li2(l-A) + Li2fl-^')+Li2fl 



'oppr 



As A — * 1, there is a manifest singularity in l4'^''^[l] since 

Ld7^(j9i,p2,P3) - log ( — ) log 



«2/ VCt4 

This double logarithm can never combine with lower point scalar integrals to form a combi- 
nation well behaved as A^^^ 0. In fact, it is easy to see from fig. ^ that the only scalar 
integrals which are available by pinching are the triangle integrals with one and two massive 
legs. These are pure poles in e and cannot be combined with the finite parts of the opposite 
box integral. There is no appropriate function which can generate the double logarithm as 
A 1 and consequently no finite function can be formed. Since the matrix elements are in 
general finite in the limit of vanishing Gram determinants, all occurences of Ldg^^ divided 
by the determinant must vanish. 

On the other hand, in D = 6 — 2e, the opposite box is not only finite as e — > as expected, 
but also as A ^ 1. This is because A^4 is effectively [A'^^^Y and its presence in the numerator 
of equation ( |2.4| ) removes the Gram determinant from the denominator. Consequently, we 
see that, 

ir ""='-''[1] = -— Ldr(Pi,P2,P3), (4.29) 

■Sl3 

which, since, 

lim Ldo(pi,p2,P3) = fsi2log( — ) +S23log( — ) -Palog f ^ 1 1 

513^0 S12S23 - P2S123 V V S12 / V S23 / \ pi J J 

+0{sl,), (4.30) 

is also finite as S13 0. 

In dealing with the all- massive triangle and adjacent box in the previous sections, we 
have constructed groups of scalar integrals that are finite both as e — > and A — > 0. For the 
opposite box, this is not easy to do since the raw scalar integral is not finite in either of these 
limits. Therefore, we follow the approach used for the two-mass triangle graph and obtain 
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the tensor integrals by direct differentiation. However, we note that the factor l/{n — D) 
associated with the differentiation formula (|2.2|) is itself divergent, implying that the 0{e) 



term in 14^^ [1] is needed - the so called e-barrier. We circumvent this, by using equation 
to construct the tensor integral with a single Feynman parameter. Explicitly, we find. 



joppr J ^ 1 [ (ai - Aa2)Ld7^(pi,p2,P3) 
(1 — A) I aiOs + Q;2a4 — aia4 — AQ;2a3 

(0204)'^ — {ci2Ci?,Y\ -^Cr / (aia4)'^A~'^ — {aia^ 




\ 04 — 03 / \ 0:4 — Aa 




(4.31) 



jopp ri ^ 1 / («2 - ai)LdgPP(pi,p2,P3) 

(1 — A) [ aiOs + a;2a;4 — aia^ — Aa2«3 



^ cr / (^2^4)" - (a2a3)' \ cr / (aia4)^A ^ - (aia3)" \ | 
\ 0:4 — aa / \ 04 — Acta / J 

This is enough to construct all the necessary integrals, since those for 02 and ai can be 
obtained by the simultaneous exchanges, ai ^ 04, 02 ^ 03 (and pi ^ ps), while integrals 
with higher powers can be obtained by direct differentiation. The only subtlety is in differ- 
entiating Ldo, where it is useful to re-express some of the logarithms as single poles in e. For 
example, 

(9Ldo(pi,j92,P3) _ (S) _ (i^) j^'^f (Qi«4)'A~^ - (a2Q4)' \ _^ 
da^ a4 — as 0:4 — Acta e \^ 0:4 y 

This helps to ensure that although the double poles in e are still divided by (1 — A), the 
single poles are finite as A — 1. In fact, when the tensor integral is multipled by a factor of e 
we also expect that the resulting single logarithms occur in groups that are finite as A — 1. 
We therefore introduce the auxiliary functions. 



Lcd„(pi,p2,P3) = -lini(exirK]) 

— — lim (a;4ail3^-'[a;4^^] — a4a2l3^'' [a^^ 



(1 - A) e-o 



(si2Lcr(pi2,P3) -p^Lc^(p2,P3)) , (4.33) 



for n = 1, . . . , 4 and. 



Lcdo(pi,i?2,P3) = 45|log i^T-^] ■ (4-34) 

^4 \P2'5l23 / 

Because these functions contain only a single power of the Gram determinant, they are not 
difficult to evaluate numerically. 
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The tensor integrals are straightforward to derive, but can be rather lengthy, for example, 



Xa 



1 - 2e (9^4 

e^ArV (^123-^12)2 {S23-Piy 

+ 7 7 - 7 2T + 7m5?^Ldo(pi,P2,P3) 

e Si3 \(Si23 - S12) (S23-P2)/ S13A 

(Sl2 - P2) 



Sl3 



-Lcdi(pi,p2,P3) - 2Lcd2(pi,p2,P3)- (4.35) 



Most of these terms are poles in e. For physical processes, the poles from different 
tensor integrals must combine in such a way that the t^l^ cancels, while the 1/e poles do 
not depend on the determinant and are ready to cancel with contributions from other loop 
configurations. In all cases, the scalar box function Ldo is associated with a single inverse 
power of . As noted earlier, for physical processes that are finite as t^l^ — >■ 0, there must 
be cancellations amongst the various tensor integrals so that no terms containing Ldo/A4^^ 
remain. We therefore choose to leave the Ldo functions exposed to facilitate the cancellation 
of these terms. 

For third and fourth rank tensor integrals, we also need to know the D = 6 — 2e box with 
one or two Feynman parameters in the numerator and the D = 8 — 2e box. The former are 
finite as e — > and are easily obtained by direct differentiation of M^l^ d-?,-2<l^ that, for 
example, 

joppD=6-2e^^^j = -^(^^i^^Ldo(pi,P2,P3)-^12Lcr(pi2,P3) 
2Si3 [ Si3 

+pM'"(P2,P3)-logf^)|, (4.36) 



S23 

while the D = % — 1e box is established using the dimension changing equation ( |2.4| ), 

yopp D=8-2e r-i 1 1 j ■512-523 — ^2-5123 y / n , /■5l23\ , / ■5l23 \ 

I4 1 = 7 \ Ldo(pi,P2,P3) - Si2log -S23log 

6Si3 \ Si3 V S12 / V S23 / 



As might be expected, both of these are finite as S13 — > and this can be seen directly from 
equation (|4.30|) . 

In summary, the situation for box integrals is very similar to that for triangle graphs. 
When the scalar integral is finite as A4 0, differentiating - or equivalently adding factors of 
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Feynman parameters - does not introduce kinematic singularities. Hence natural groupings 
of box and triangle integrals arise that are finite as A4 —>■ 0. Furthermore, the infrared 
and ultraviolet singularities can be isolated easily. Although we have explicitly worked 
through a subset of kinematic configurations relevant to certain QCD processes, this method 
is systematic and can be applied to processes with more general kinematics (and particularly 
non-zero internal masses). 



5 Five point integrals 



In this section we consider five point integrals with only one external mass. The outflowing 
lightlike momenta are denoted pi, i = 1, . . . ,4 while the fifth = — P1234 is massive, p^ ^ 
as shown in fig. |^. The auxiliary momentum is then, 

= -(1 - xi)p'^ - (1 - a;i - 0:2)^2 - (3^4 + 3:5)^3 - ^5P4- 
We can make the choice, 

ttlttS 51234 = -A, 



aia3Si2 
aia4Si23 



tt2tt4 S23 

a2a5 S234 
a3«5 S34 -- 



with, 



^1 - A) = TT ai(si23S234 - S23S1234) 

"3 ,ti 



As in the opposite box integral, A is an extra kinematic variable that ensures the at are 
independent. It is the same variable that occurs in the third pinching which forms an 
opposite box configuration. The matrix rjij is given by, 

/I -1 1-A 
-1 1 A-1 
[l-Xf 
A-1 
1-A 



Vij 



1-A A 

1 1 - 
-1 



- 1 
2A 



1 - 2A 
A-1 
1 

-1 



-1 \ 
1 

1-A 
-1 
1 



and the normalisation factor is. 



A. 



The 7i are rather lengthy, but can be read off from rjij. The scalar pentagon integral is by 
now well-known in D = A 0, ^ and in D = 4 — 2e @, |^, and can be written in terms 
of these variables as, 

h[i] = -^j:iJr\i]+o{e). (5.1) 



^ 1 
3 m=l 
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p^+P2+P3+P4 







PI+P2 




P3+P4 




Figure 7: The pentagon graph and each of the five pinchings obtained by omitting the 
internal fine associated with am for m = 1,2,3,4 and 5. For box pinchings, the entering 
momentum is fixed by momentum conservation. 



The five pinchings and the momenta associated with each is illustrated in fig. |^. The limit 
A^5 — *■ corresponds to the vanishing of the Gram determinant associated with the m = 3 
pinching. The scalar integral for this pinching is not well behaved in this limit and so should 
not be expected to combine with the other pinchings. Therefore, we separate l5[l] according 
to the pole structure in e and N^. We identify the function Lei which is finite as both ^ 
and e — >• and does not depend on the opposite-box pinching, m = 3, plus a combination of 
scalar box integrals containing all the infrared poles and the remaining ^ singularities. 



l5[l] 



^ a^a.K^ Lei - ^Tslf [1] - a2a4 (if [1] + if [1] 
-^as (ailf [1] - a2lf [1] - «4lf [1] + a,lf\l]) , 



where = (1, —1, 0, 1, —1) and. 



Lei 



2m 



5 

E 

m=l 



Kir 



(5.2) 



(5.3) 



Applying equation (|2.7| ) and noting that If ^ is both infrared and ultraviolet finite, the 
integrals with one insertion are also determined in terms of box integrals in Z) = 4 — 2e 
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dimensions, 



h[a^] = -^j:V^J^\l]+0{e). (5.4) 

^^^5 rn=l 



Making the same separation as before, 

1 ' 



h[xi] = [asaiKi) Lei - tttt H "i"™ {Vim - i^ii^m) + C'(e), (5-5) 



where the bracket multiplying the scalar box integral is always proportional to N^. 



Similarly, using the usual formula for two Feynman parameters ( p.8|) and concentrating 
all of the As dependence into Cij, where, 

= a.ajC^j = ^ \ r]ij - ^ I > (5-6) 



we find. 



1 5 1 5 



4^^5 m=l ^^^5 m=\ 

Asirni] 



5 



— E if =^'('") [1] - r^.Jr^ [a.] I . (5.7) 



To simplify the non-Qj- terms (and eliminate one power of A^^s), we have rewritten some of 
the variables appropriate to the pentagon integral in terms of those appropriate to the box 
pinchings, r/lj"^ and 7j-'"'', using the relations 

(m) '^mm'Ji ^im'ym (m) VmmVij VimVjm 

~ 2iV^ ' ~ 2iV5 ■ 



Bern, Dixon and Kosower have shown P, |^ that I5 drops out of the calculation of any 
^auge theory amplitudes by using the identity H], 



^ 1 



E grgfQ+i,,+i = -^/|:f , (5.8) 
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where g^^-l^^ represents the metric tensor in D = 4. Since, 

= j: l5[x.+ix,+i]grgf - ht'lg'-''], (5.9) 

a// terms proportional to Cjj in l5[ajaj] can be moved into the existing g^^^^ piece. Inspection 
of equation ( |5.7| ) indicates that the correct term to reshuffle is — A5i5[l]/A'"5 rather than 
merely Retaining only the Cjj terms, 

terms = -^I^='^[l]^^'^^ - ^^^^"'^^^ ^ gf^f Q+U+i 

= 4^E7™ir''^™^[lK''^^+0(e). (5.10) 

5 m=l 

As expected the terms precisely cancel. Here the finiteness of and If ^'^[l] has 

been used to ensure that this term generates only 0{e) corrections when replacing g'^l^^ with 
the dimensionally regularised g^^^^^K The remaining piece should not contain any kinematic 
singularities associated with lower-point Gram determinants, in particular N^, as it originates 
from the well-behaved if This is indeed the case and we write, 

^^2 = ^j:iX=''^"M- (5.11) 

The pentagon integrals with three insertions can be obtained by direct differentiation of 
equation (|5.7| ). We find, 

D=6,(m) 
4 Ffc. 



5 r7i=l \ \ ^mm / 



(m) (m) 

VimVjk ~ Vij Vkm \ ^D=6,{m)i^-i f(m) 



I4 '^""-"[1] - r]jml4 [o-iak] + cyclic i,j,k.(5.12) 



Since the Cjj term is obtained by differentiating Cjjif '^[l], it must be finite as both ^ 
and e ^ and we introduce the finite function, 

Lesfc = — - V (afcam77fcmlf "^'^""^[1] + 27mlf "^''•"''[xfc]) . (5.13) 



36 



Recalling that, 



5 

= - H hlxi+iXj+iXk+ijqtqj'qk" 

i,j,k=l 

+ 1 {E ^t'[^k+i]9'''^^q^k' + cyclic i,j,k| , (5.14) 
and keeping only the Cij terms in equation (|5.12|) we find, 

g,.Hq^^^ terms = \lt'[xk+i]g'^'^qi:' - {it'l^k+i] + Le^k) E c.+i,,+igrgf 



= -lLe3,r/'^'^^g^'= + 0(e). (5.15) 

As in the previous case, the finiteness of the coefficient of g'^^^^ has been used to promote it 
to the full D = A — 2e metric tensor. 

The tensor integrals with four and five insertions may be obtained by further differentia- 
tion, and the same trick used to rewrite the Cjj terms as a contribution to the metric tensor 
structure. In this way, all vestiges of the pentagon in D = 6 — 2e and higher dimensions can 
be removed, along with the inverse powers of A5. 

The functions introduced in this section. Lei, Le2 and Le3j, only contain a single power of 
the Gram determinant. Consequently, they can be evaluated without numerical problems. 
To illustrate this, we choose a representative phase space point, S1234 = 1, S123 = 0.4, S234 = 
0.3, Si3 = 0.1 and use the variable S23 to control A. The sixth variable S24 is chosen to lie 
within the physical region defined by the pentagon Gram determinant, A5 < 0. Fig. |] shows 
the functions Lei and Le2 together with Le35. In each case, we see that the A — >■ 1 limit is 
smoothly approached indicating that the function is intrinsically well behaved in that limit. 



6 Conclusions 



In this paper we have developed a new strategy for evaluating one-loop tensor integrals. 
It avoids the usual problems associated with the presence of Gram determinants. Such 
Gram determinants arise when the tensor integrals are expressed in terms of the physical 
momenta and generate false singularities at the edges of phase space. In addition to creating 
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10 - 



5 - 




10-'^ 10-" lO"'" 10"^ 10" 10'' 10""^ 10"'' 10'* 10'^ 10'^ 10"' 10° 

1-1 

Figure 8: The finite functions for the pentagon graph with one external mass and evaluated 
in double precision Fortran as a function of 1 — A. Because the functions only contain single 
inverse powers of (1 — A), no numerical problems are evident. 



numerical instabilities, they tend to increase the size of the one-loop matrix elements. Our 
approach is to construct groups of scalar integrals which are well behaved in the limit of 
vanishing Gram determinant (A„ — >■ 0), and which can be evaluated with arbitrary precision 
by making a Taylor series expansion in A„. In fact such combinations arise naturally by 
either differentiating with respect to the external parameters - essentially yielding scalar 
integrals with Feynman parameters in the numerator - or by developing the scalar integral 
in D = 6 — 2e or higher dimensions. Evaluating these new integrals is straightforward - they 
are just linear combinations of the known scalar integrals in Z) = 4 or D = 4 — 2e. As such, 
they combine the dilogarithms, logarithms and constants from different scalar integrals in 
an extremely non-trivial way. As a bonus other spurious kinematic singularities are also 
controlled - they appear in the denominator of the finite functions, which are well behaved 
in the singular limit. Although the number of basic functions has increased, the number of 
dilogarithm evaluations has not, since the functions are generated recursively. Furthermore, 
because the Gram determinant singularities are not genuine, by grouping integrals in this 
way, the expressions for one loop integrals are compactified. 

To illustrate the approach for specific integrals, we have applied the method to 3-, 4- and 
5-point integrals where the internal masses have been set equal to zero. These tensor integrals 
are relevant for a range of QCD processes where the quark and gluon masses are negligible. 
For more general processes with arbitrary internal masses and external kinematics, the rel- 
evant combinations of scalar integrals can be obtained using equations ( |2.4| , |2.7| - |2.10| ). As a 
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by-product we have shown how all the Gram determinants associated with pentagon graphs 
can be eliminated. 
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A Triangle integrals 



In this appendix, we collect together explicit forms for the triangle graphs which appear 
as building blocks in the box graphs. We fix the kinematics according to fig. 1, so that 
momenta pi and p2 are exiting, with determined by momentum conservation. There are 
two distinct cases, according to whether one or two of the external masses are zero. 



A.l The one-mass triangle 

Here we provide explicit results for the case pf = P2 = ^- There is an additional symmetry 
under the exchange pi <-> p2 and Xi <-> X3. Insertions involving X2 can be eliminated using 

Xi+X2 + Xs= 1. 

The D — 4: — 2€ scalar and tensor integrals are given by. 



I-[l] = ^i^^, (A.1) 



S12 

ll"^[x,]=ll"^[xs]=2ll"^[xl]=2ll"^[xl] = -J-(t^ + 2]cr, (A.2) 

I'rix.xs] = 3ll"^[xlx,] = 3ll"^[x,xl] = (A.3) 

^Si2 

rH = ir[xil = -3^^ (^^ + f ) ^r. (A.4) 
while, the necessary integrals in D = 6 — 2e dimensions read, 

I3 [1] = -2 1 1 + 3lcr, (A.5) 

I3 [xil^h N = -g I +3)^^- (A.6) 



A.2 The two-mass triangle 

When two of the external legs are massive, pi — but S12, pi 7^ 0, the divergent integrals 
read. 
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If 



-2Lcrbi,p2)- 
4 



(Ai 



4 
Pi 



2(si2-p?) 



cr- 



e 2(512 



(A.9) 



y^'^ ^^^'^^^ + 2(.i2-1>?)^"(.12-pD^'^ 



.-^12 



+ 



4 
Pi 



3(S12-P?) 2(S12-P?)2 (Sl2-P?)3 



while, 



If [xi] = 2lf [xf] = 3lf [xf] = Lcr(pi,P2), 



-Lc2™(pi,p2). 



The functions Lcf (pi,p2) are defined by equations ( p.l9| , p.20D 
The D = Q — 2e dimension integrals read, 

1 



2m,D=6-2e 



I 



j2m,_D=6-2e 



[1] = - pfLcr(Pi,P2) 



I 



2m,_D=6-2e 



[xi] = - (^pfLcr(Pi,P2) 

[X3] = ^ (p?Lcf (pi,P2) 



e 

^-^12)' 
e 

;-si2)^ 



3 cr 



cr 



3 



cr, (A.IO) 

(A.ll) 
(A.12) 

(A.13) 



(A.14) 
(A.15) 
(A.16) 



The corresponding integrals for the case p\ = 0, ^3 7^ can be obtained by using the 
above formulae with the substitutions. 



Pi P2, 



Xl ^ X3. 



(A.17) 



The limit — > may also be safely taken. For example, using the fact that in this limit 
Lcf (pi, P2) ^1 / S12 and observing that all of the triangle loop integral contributions 
proportional to trivially drop out in eqs. (A.8)-(A.10), we see that. 
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B Box integrals 



Here we collect together explicit forms for the some of the finite functions appearing in the 
box integrals. We fix the kinematics according to fig. 2, so that momenta pi, p2 and are 
exiting, with p4 determined by momentum conservation. 



B.l The adjacent two- mass box 

We first focus on the case where two of the adjacent external legs are massive, pi 7^ and 
pj^O. 

For integrals with a single Feynman parameter in the numerator, we have, 

Ldi(pi,P2,P3) = - (Ldi5(pi,J92,P3) + Lco(pi,P23)) , (B.l) 

where the box function in D = 6, Ldis is given by equation ( [4.3|) . When there are two 
Feynman parameters in the numerator, we eliminate X3 using xi + X2 + X3 + X4 = 1 so there 
are only three relevant functions. Explicitly, we find, 

2 

Ld2l(pi,P2,P3) = (3Ld25(Pl,P2,P3) + Lcis(pi,p23)) - Lci(p23,Pi), (B.2) 

S12 

Ld22{pi,P2,P3) = '^^^-^ ^ (3Ld25(Pl,P2,P3) + Lci5(j9i,J923)) 

■Sl2"S23 

-^Lc^(pi,p2) + ^Lci(p23,Pi), (B.3) 
S23 S23 

Ld24(Pl,P2,P3) = ^ (3Ld2s(pi,P2,P3) + Lci5(pi,p23)) 

■Sl2'523 

_£11LC^(J912,P3) + ^Lci(j923,Pl). (B.4) 
S23 -523 

The all massive triangle integral function in D = 6 — 2e, Lcis is given by equation ( |3.11| ), the 
box integral function in D = 8 — 2e is given in equation (|4.7| ) while the remaining triangle 
functions are given in Appendix A and equations (^.4|-p.9D. The functions for adjacent box 
integrals with three insertions of Feynman parameters contain the box in D = 10 — 2e ([4.8|) 
and triangle in Z) = 8 — 2e ( p.l4| ). All integrals can be obtained in terms of the following 
four functions, 

/ 12 

Ld3ii(pi,P2,P3) = - — (5Ld35(pi,p2,P3) + LC25(Pl,P23)) 

\^12 

+ ^^i^^Lc3(p23,Pl) + ^LC2(P1,;^23)V (B.5) 
S12 Si2 J 
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\2(p'^ -512) 

S12S23 

H ^Lc3(j923,Pl) - — ^LC2(pi,P23) 

S12S23 S12 

^^'Lc^(pi2,P3), (B.6) 



2S23 

Ld322(pi,P2,P3) = -( ^^^^^T-2-^^ (5Ld35(pi,P2,P3) + LC25(Pl,P23)) 

V ■512^23 

Si23(Sl2Sl23 +P?Si23 - 2p?Si2) . ^ S12 2m/ N 

^ 2 LC3(p23,Pl) - 7T— LC2 (Pl,P2) 



S12S23 2S23 

, -512(512 — Si23)y 2m/ N , -5123(5123 — 2512)^ , . 

+ TT^ (pi2,P3) + LC2 Pi,P23 

2S23 -S12-S23 



(P?Sl2 + -S12-S23 -Sl2-Sl23j y 2m/ 



+ (pi,P2) , (B.7) 

'23 / 



(12(s — n^)^ 
^ ^ (5Ld35(pi,p2,P3) + LC25(Pl,P23)) 
-S12-S23 

,Pl(Pl-2Sl2) . Pi(Pl-^12) y / X 

H LC2(pi,P23) H 2 Lc3(p23,Pl) 

-§12-523 -512-S23 

-^i^%^Lcr(pi2,P3) + ;^Lcr(pi2,P3)V (B.8) 



2S^3 ^ ^^^^'^-^^ 2S23 



B.2 The one-mass box 



For this kinematic configuration, only and there is a 'flip' symmetry, so that functions 
related to the parameter X4 are obtained from those related to xi, by pi < — > p^. The box 



integrals in higher dimension are given by equations (|4. 181 - ^^211) , and, 

Ld^(pi,P2,P3) = -Ld}^(pi,P2,P3). 

For two and three insertions, we flnd. 



and, 



Ld2r(pi,P2,P3) 

L4r(Pl,P2,P3) 

Ld3n(pi,P2,P3) = 



1 



(6Ld^^(pi,P2,P3) + S23Lcr (P23,Pl 

-S12 ^ 

Ld}™(pi,P2,P3)-Ldlr(pi,p2,P3), 



.^f60Ld^^(pi,P2,P3) + ^S^3Lc^(p23,Pl) 
-^12 \ ^ 



(B.9) 

(B.IO) 
(B.ll) 
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+ ^S12S23LC^(P23,P1)), (B.12) 



S12S23 \ ^ 



+ ^4Lcf"(pi2,P3)j, (B.13) 
L4^2(Pi,P2,P3) = L4ri(Pi,P2,P3) + Ld^(pi,p2,P3)-2Ld^r(pi,p2,P3)- (B.14) 



44 



C Limits 



In this section, we collect together suitable expansions of the various functions presented 
in this paper in the hmit A„ — > 0. These expressions represent the leading term in the 
expansion of the functions as a Taylor series in A„ and should be evaluated for A„ < S, 
where S can be determined numerically. Typically, 6 ~ 10~^A™'^^ where A™*^^ is the largest 
value the Gram determinant can achieve. In general, for a given numerical precision, acc the 
numerical problems occur when A ~ {acc)^/^ where N is the number of Gram determinants 
in the denominator of the function. 



C.l The three-mass triangle 



In the limit that A3 0, we have. 



1 



LCo(pi,P2) Tnr2 Pl(«12+P2-Pl)l0g^ +P2(«12+P?-P2)logl 2 



(C.l) 

Lci5(pi,P2) ^ UpIpIsu (^^^^'' ^ ~ ( ) ^ ^ ~ ( 



+2p?p2^,2j, (C.2) 
+^^{su+pI+pI)), (C.3) 



1680^ ( P?(^i2 + - pI) log ( + P^(.i2 +pI- pI) log 



+^^if^(^?2 +P\+P'. + Si2pI + ^i2P^ + pIpI?! , (C.4) 

Lci(pi,p2) - 3pf(.,2+pi-pf) (^) +2p^-ptLco(pi,P2)) , (C.5) 

LC2(P1,P2) - 6pl(.,2 +p! -Pi) (~^^^^'^^ +Pl -P2)LC3(P1,P2) 

+Sl2 -P2 + 3PiP2LCl(pi,P2) -PiLCi(p2,Pl) +P2l0g (^j j' i^'^) 
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LcM,P2) - 5pi{s,,lpl-pl) ((^^^-^^)^"g ) +Pi-2ptLc,(p„p,)j , (C.7) 

+S12 -P2 + 5p?P2LC3(Pl,P2) - 4p^Lc2(pi,P2) +P2log ) ' '^^'^^ 

LC5(P1,P2) - 7pf(.i2+pi-pf) (^'" " ) ^ " 3ptLC3(pi,P2)) . (C.9) 



C.2 The adjacent two-mass box 

In the limit A4 ^ 0, we have, 

Ldo(pi,P2,P3) "2 ( «123 +Pl - S23 ^-^j Lco(pi,P23), (C.IO) 

Ldi5(pi,P2,P3) ^ ( +Pl - S23 - ) Lci5(pi,P23) 

•S12S23 \\ S\2 J 

2 



S23 / \Si2 J Z \pl 

2 2pfsi23^ 



Ld2s(pi,P2,P3) ^ 77^ -5123+^1 ^ S23 LC2s(pi,P23) 

IZS12S23 \ \ •S12 / 

•^23 i„„ / -5123 \ 



--log 

2 V S23 / 



+ .Ll0g(^)-|l0gP? +.12.23 , 

S23 / V S12 / 2 V Pi / / 

(C.12) 



Ld35(Pl,P2,P3) 180.12523 ( (^^^^ ^ ~ ~ ~^^) ^^^^^^^'^23) 

.3 /„ X /„ X ^6 /„ ' 

.123 



, ^23 1 /Sl23\ , 3 1 /Sl23\ Pi , 
+ ^log +Sl2l0g -—log, 2 

2 V .23 / V .12 / 2 \pi 



+ ^.12.23(Sl23 + S12 + .23 + Pl)j , (C.13) 



Ld4s(Pl,P2,P3) 840.12.23 ((^^^^ "'"^1 ~ ~ ~^^) ^^'^^'^^l'^23) 



, ^3 , /Si23\ , 4, /Sl23\ Pf, ( Sl23 

1 



.23 / V .12 / 2 \ pi 



+ 2*12.23(^123 + *12 + *23 +Pl + '5l23'5l2 + .123*23 
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+ s,,pI + s,sPi + ^)) . (C.14) 

In this last equation, we have used the finite part of the three mass triangle graph in D = 
12 -2e. Forp?, pi ^ 0, 



+Pl{si2+pl-pl)l0g 



Pi 



+ ^-^^{pt + pt + sl, + pIs,, + pIs,, + pIpI)] \ . (C.15) 



C.3 The one- mass box 



Finally, in the limit — > 0, we have, 

'Ld]r{pi.P2.Pz) ^ 0, (C.16) 

Ldig'(Pl,P2,P3) ^ ^(S23l0g(^)+Si2l0g(^)), (C.17) 

•812523 V V *23 / V ^12 / / 

Ldg(j)l,P2,P3) ^ B;^(*23l0g(^)+i!;2l0g(^)+Sl2i>23), (C.18) 

+ ^S12S23(S123 + Sl2 + ^23)^ , (C.19) 

1 / 



S12S23 V V S23 / V S12 / 

I 1 (2,2,2, I I -812-523 \ 
+ g'8l2'523 I ^123 + -812 + *23 + '5l23'Sl2 + 5123*23 "1 ^ I 



(C.20) 
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